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Abstract 

With its ability to look at bright galactic X-ray sources with sub-millisecond time 
resolution, the Rossi X-ray Timing Explorer (RXTE) discovered that the X-ray 
emission from accreting compact stars shows quasi-periodic oscillations on the dy- 
namical timescales of the strong field region. RXTE showed also that waveform 
fitting of the oscillations resulting from hot spots at the surface of rapidly rotating 
neutron stars constrain their masses and radii. These two breakthroughs suddenly 
opened up a new window on fundamental physics, by providing new insights on 
strong gravity and dense matter. Building upon the RXTE legacy, in the Cosmic 
Vision exercise, testing General Relativity in the strong field limit and constraining 
the equation of state of dense matter were recognized recently as key goals to be 
pursued in the ESA science program for the years 2015-2025. This in turn identi- 
fied the need for a large (10 m 2 class) aperture X-ray observatory. In recognition 
of this need, the XEUS mission concept which has evolved into a single launch L2 
formation flying mission will have a fast timing instrument in the focal plane. In 
this paper I will outline the unique science that will be addressed with fast X-ray 
timing on XEUS. 
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1 Fast time variability and strong gravity 



The X-ray emission from accreting compact stars (neutron stars and black 
holes) has been shown to vary on (sub)-millisecond timescales, i.e. on timescales 
comparable to the dynamical timescales of the innermos t regions of the ac- 
cretion disk. This variability was discovered with RXTE (|Bradt et all Il993j) 
just ten years ago in the form of kilo- Hz quasi-periodic oscillations (QPOs) in 
Fourier power density spectra (see the review by van der Klis (2005) and Fig. 
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Fig. 1. Power density spectrum of Sco X-l in a vFv-like representation, showing 
two kHz quasi-periodic oscillations (data courtesy of Michiel van der Klis). 

1). There is a wide consensus that these QPOs probe the motion of matter un- 
der strong gravity with the signal originating from within a few Schwarzschild 
radii of the compact star. This is a region where the spacetime curvature is 
extreme, orders of magnitude larger than that sampled in weak field tests (e.g. 
Gravity-Probe B). This is a regime where the most dramatic effects of GR are 
expected, and where fundamental predictions of GR, such as the existence of 
an innermost stable circular orbit (ISCO) have yet be tested. With RXTE, 
X-ray timing has thus become a complementary tool to X-ray spectroscopy, 
polarimetry and gravitational wave measurements to probe strong field GR. 
It has even been claimed that studyi ng black hole variability might be used 
to test alternative theories of gravity ( Psaltisl . 2004 ). 



1.1 Overall QPO properties 



It is is beyond the scope of this paper to review the properties of kHz QPOs, 
and the reader is referred to van der Klis (2005) for a recent review and X- 
ray Timing 2003 - Rossi and Beyond by Kaaret, Lamb and Swank (2004, 
AIP Conference Proceedings, Vol. 714, Melville, NY: American Institute of 
Physics) on which the present paper draws heavily on. Here I will focus on 
their general properties, relevant to some of the points discussed below. So 
far, kHz QPOs have been detected from a bout 20 neutr on star low-mass X- 
ray binaries (about 10% of known systems) ((Swank . 2004 ). and from a handful 
of binaries containing a black hole candidate. In general two kHz QPOs are 
observed. The properties of the QPO (frequency, amplitude, width) vary in a 
complicated way with parameters such as the source luminosity. For neutron 
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stars, the QPO frequencies can vary by at least two orders of magnitude with 
the amplitude decreasing with increasing frequency. There is a trend for the 
QPO fractional amplitude to increase with energy, reaching in some cases 
20% at ~ 20 keV. In black holes, QPOs are detected closer to the sensitivity 
limit of RXTE with typical amplitudes of 1%: they are transient, show stable 
freque ncies and are de tected in the hardest X-ray bands (typically above 10 
keV) (|R.emillardl . l200l . 



1.2 QPO interpretation 



The true nature of the underlying signal is unknown and there is a wide 
range of possible models to be considered, mostly focussing on predicting 
QPO frequencies. Because the disk is a natural source of periodicities, from 
orbital to epicyclic motions and disk oscillations, it is generally agreed that 
QPOs originate in the cool (~ 1 — 2 keV) accretion disk. This implicitly 
assumes that the oscillator must survive the strong damping forces thought 
to be present in the disk and that an amplification mechanism operates to 
account for the fact that QPOs show larger amplitude at higher energies (~ 
10 — 20 keV). For neutron star QPOs, the main constraint for the models 
so far comes the observation that the frequency separation between the twin 
QPOs is close to the spin frequency or half the spin frequency of the neutron 
star (where that is known). This strongly suggests that the spin is involved 
in the generation of the QPOs. For black holes, the main constraint for the 
models comes from the observations that QPOs, when observed in pairs, have 
commensurate fr equencies, with ratios close to small integer ratios (e.g. 3:2) 
( Remillardl . 2004 ) , lending support to the idea that a resonance mechanism is 
involved in the make up of the QPOs. 



Some models (not all of them consistent with the constraints above) relate 
QPO frequencies directly characteristic frequencies of test particles moving 
in the stron gly curved spacetime (General R elativistic orbital and epicyclic 
frequencies ( Abramowicz and Kfuzniakl . 2004) ) , others link the QPO frequen- 
cies to the fr equency of orbitin g clumps interacting with radiation from the 
neutron sta r dMiller et ah . 1998j ). Some models involve global disk oscillations, 
relativistic (|KatoL l200ll T or not ([Titarchuk and Woodl 120041 ) . Yet others as- 
sociate QPOs with strong field GR e ffects, e.g. relativist ic dragging of inertial 
frames ( Stella and Vietrif Il999f ). see ( van der Klisl . 2005| ) for a more complete 
list of references. 



There is however a trend in the models to associate QPOs with resonance 
phenomena, involving General Relativistic frequencies. As an example, in the 
accreting millisecond pulsar, SAXJ1808-3658, it has been proposed that the 
two QPOs (whose frequency difference is just equal to half the neutron star 
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spin frequency) are generated at a radius in the disk where the difference be- 
tween the general relativistic orbital frequency and radial e picyclic frequency 
is equal to half the spin frequency ( Wiinands et al. . 2003j ). Similarly, to ac- 
count for the 3:2 frequency ratios observed in black hole systems, a parametric 
resonance concept has been put forward, in which the QPOs are produced at 
a radius in the disk where two of the three general relativistic frequencies (or- 
bital, vertical and radial ep icyclic) have commensurate value s, matching the 
observed QPO frequencies (jAbramowicz and Kluzniakl . 12004). Although this 
is still very much under discussion, it is interesting to note that the latter two 
models link QPOs directly to general relativistic frequencies. 



1.3 The innermost stable circular orbit 



The ISCO is one of the key predictions of strong gravity GR stating that 
there exists a region around sufficiently compact stars within which no stable 
circular orbital motion is possible. In a Schwarzschild geometry, the radius of 
the ISCO is 6GM/c 2 , which is larger than the radii of neutron stars, deduced 
from models c onstru cted with most modern high-density equations of state 
f Akmal et all ^998). Interestingly enough, claims have already been made 
that the QPO properties have revealed signatures of the ISCO. Prior to the 
launch of RXTE, it was suggested t hat the ISCO might i nduce a frequency 
cutoff in the power density spectrum ( Kluzniak et al. , 1990j ) . After the discov- 
ery of kHz QPOs, it was proposed that signatures of the ISCO could include 
a frequency saturation of kHz QPOs with increasing mass accretion rate, or a 
drop in the qual ity factor and ampl itude of kHz QPOs as they approach a lim- 
iting freque ncy dMiller et al. , 1998| ). A saturation of the QPO fr equency with 
count rate ( Zhang et al. . 1998| ) or inferred mass accretion rate ( Bloser et al 



2000[ ). as well as a sudden drop of coherence and amplitude at some critical 
frequency ( Barret et all l2005bf ) have now been reported. Although the inter- 
pretation of these results is still debated, this demonstrates the potential of 
fast X-ray timing for testing the existence of the ISCO. Note that identify- 
ing securely a frequency with an orbital frequency at the ISCO would yield 
directly the mass of the neutron star. 



1.4 QPO and the black hole spin and mass 



Depending on the interpretation, parameters such as the spin and the mass of 
black holes can be constrained directly from measuring the QPO frequen- 
cies. One example is given in Fig. 2 where the most model independent 
constraint is used. It simply assumes that the QPO frequency is smaller 
than the orbital frequency at the ISCO (the latter being a function of the 
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Fig. 2. Left) The azimuthal (/<a), radial epicyclic (k), and vertical frequen- 
cies at different radii around a black hole with spin parameter a/M = 0.2. Right) 
Dashed Lines: the dynamical measurement of the black-hole mass in the source 
GRO 1655—40. Solid Line: the minimum spin parameter a/M for each black-hole 
mass, for which the observed 450 Hz QPO can be produced as a Keplerian frequency 
of a stable orbit. According to this argumen t, the b l ack h ole in GRO 1655—40 is 
spinning with a/M > 0.2 (Both figures from rtPsaltislEooi l) 



black hole mass and spin). More stringent constraints on the black hole mass 
and spin can be o btained from models involving resonance betw een rela- 
tivistic frequencies (|Psaltisl . l2004t Abramowicz and Kluzniak . 2004 ). All con- 
straints point to maximally spinning black holes. As an application, because 
of the 1/M dependency of all general relativistic frequencies, it has been pro- 
posed that measuring QPOs from Ultra-Luminous X-ray sources, whose na- 
ture is a mystery, might help in c onstraining the mass of the compact object 
(|Abramowicz and Kluzniakl . 120041 ) in these systems. 



2 Fast X-ray timing and matter at supra-nuclear density 



Besides providing a laboratory to test GR, neutron stars have also other in- 
teresting properties; some may be powerful emitters of gravitational waves at 
their millisecond spin periods, some have an extreme magnetic field in which 
strong-field QED effects can be studied. But the overall reason why neutron 
stars have stimulated much work is because the physical state of the matter 
which makes up their cores is still a mystery. The density of a neutron star can 
reach up to 10 times the density of a nucleus; besides neutrons and protons, 
at such high densities the physics of strong interactions predicts the existence 
of hyperons - particles that contain strange quarks - pions and kaon conden- 
sates, or even free quarks. Neutron stars offer the only possibility of detecting 
such exotic matter in the bulk, and resolving the question of the most stable 
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form of matter (baryons or free quarks). Since the structure of neutron stars 
is set by the, yet unknown, equation of state of matter at these high densities, 
measurements of neutron stars masses and radii can provide information not 
only on the physical makeup of neutron stars, but also on the nature of the 
interactions between the particles that constitute them. The equation of state 
of nuclear matter is at the heart of many fundamental problems in physics, 
from nuclear physics to condensed matter physics, including superconductivity 
and super-fluidity. Fast X-ray timing and high resolution spectroscopy (e.g. 
redshifted lines emitted at the neutron star surface) holds great potential for 
constraining the physics of strong interactions. 



2. 1 QPOs and the equation of state of dense matter 



If QPOs arise from stable orbital motion (the highest frequency QPO is gener- 
ally associated with an orbital frequency in most models), the frequency yields 
the orbital radius at which it is produced. Hence this orbital radius limits the 
mass and the radius of the neutron star ( Miller et all fl998h . This is because 
two types of constraints can be set. The first one uses the fact that no QPOs 
are expected to be produced within the ISCO, hence the orbital radius must 
be larger than the radius of the ISCO (6GM/c 2 in a Schwarzschild geometry). 
The second constraint simply states that the radius of the neutron star is 
smaller than the orbital radius. So far, the highest QPO frequency is 1330 Hz, 
limiting the mass and radius of the neutron star in the bottom left area of Fig 
3. 



2.2 Waveform fitting of neutron star X-ray oscillations and dense matter 



The discovery of the first accreting millisecond pulsar (a handful are known 
today) and transient oscillations in type I X-ray bursts in a dozen sources 
have provided us with additional tools to constrain dense matter. This is 
because these pulsat ions ar e likely caused by hot spo ts rotating at the surface 
of the neutron star ( Munol . 2004 ; Strohmaver . 2004 ). The emission from the 
spot will travel in a strongly curved space time before reaching the observer, 
and will be subject to a variety of effects, including relati vistic aberration, 
Doppler boosting and strong field gravitational light bending ( Poutanenl . 2004 : 



Strohmaverl . 120041 ) . The more compact the star, the smoother will be the pulse 



profile. Fitting of the pulse profile of the oscillations with generic light curves 
computed under specific assumptions (spacetime metric, geometry, emission 
pattern and spectrum of the spot,....) has already demonstrated the potential 
of this approach to constrain the compactness of the st ar, especially when 
the harmonic content of the oscillations is also detected ( Strohmaverl . 120041 : 
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Fig. 3. left) Constraints from orbital frequencies. The 1330 Hz curve is for the 
highest kilohertz quasi-periodic oscillation frequency yet measured. This curve is for 
a non-rotating star; the constraint wedge would be enlarged slightly for a rotating 
star. The solid lines are mass-radius curves for different representative high-density 
equations of state. Right) Constraints on mass and radius for hypothetical detections 
of a 1500 Hz QPO and a 1800 Hz QPO, if they are identified with an orbital 
frequency. At 1500 Hz one expects signatures of the ISCO to be present; a detection 
of 1800 Hz would present strong difficulties for standard nucleonic equations of state. 
Both figures taken from Miller (2004). 

Bhattacharwa et al. . 20051 ). 



3 From RXTE to XEUS: From 0.7 m 2 to ~ 10 m 2 



As stressed above, we are at the beginning of a new era, which means that the 
full potential of the RXTE discoveries is difficult to anticipate. It is however 
clear that fast X-ray timing holds promises to address some of the most im- 
portant issues of modern astrophysics, such as strong gravity, accretion disk 
physics, and the physics of the strong interactions. The obvious path forward 
will therefore require an instrument with a significant increase in sensitivity 
compared to RXTE. For timing, the significance of a QPO detection scales 
proportionally with the number of photons collected per unit time, hence the 
with collecting area of the instrumentation. Given that RXTE has a collecting 
area of ~ 0.7 m 2 , a ~ 10 m 2 class X-ray telescope (with say 3 m 2 effective 
area at 6 keV) has the potential to provide an improvement in sensitivity for 
timing measurements by more than one order of magnitude over RXTE. Such 
an instrumentation would allow breakthrough science in several ways. The 
most important ones are listed below. 
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3. 1 Breakthrough 



It is already clear that with RXTE, we have only seen the tip of the iceberg of 
the high frequency phenomena around compact stars. A factor of 3 to 4 weaker 
signals in terms of RMS amplitude will be detected with a 10 m 2 instrumenta- 
tion. Side-band frequencies, predicted in some models, will be detected, whose 
inter-dependence will help in revealing their nature. Measuring simultaneously 
three frequencies could establish their identification with GR frequencies, if for 
instance they vary in the way predicted. New analysis techniques, such as the 
one based on the bicoherence, which measures the degree of phase coupling 
between oscillations at diffe rent frequencies, might be sensitive enough to dis- 
criminate between models (M accarone and Schnittmanl . l2005h . Furthermore, 



because the RMS amplitude of the QPO decreases with increasing frequen- 
cies, QPOs will be detected at frequencies higher than currently achieved with 
RXTE. This would allow a search for the frequency ceiling expected at the 
ISCO and, provided that the orbital nature of the signal is demonstrated, 
would restrict critically the region of acceptance for neutron star masses and 
radii, hence allowing some proposed equations of state to be eliminated (see 
Fig 3). 

A second breakthrough would be the determination of the underlying modu- 
lation mechanism of the signal. In general, it is simply assumed that a QPO 
signal is made of sine-wave damped exponential shots, with a time constant 
(also called coherence time) which relates to the measured width of the QPO 
in the Fourier spe ctrum. Coherence tim e of the order of 0.1 second have been 



inferred this way ((Barret et all l2005allbl ) . Because this estimate is computed 
by averaging the signal over thousands of cycles it should be considered as a 
lower limit (frequency drifts will artificially broaden the signal, thus reducing 
the inferred coherence time). A 10 m 2 class instrumentation will enable to 
study kilo-Hz QPOs on timescales comparable to their coherence times. This 
will allow one to test whether the signal is persistent or transient, intrinsically 
aperiodic or made of a superposition of more coherent pulses. The true lifetime 
of the underlyi ng oscillator may wel l be one of the most stringent constraints 
to the models ( Barret et al. . 2005bh . Depending on the nature of the under- 



lying signal, time-resolved spectroscopy over its coherence time may become 
possible, yielding information on the radius where the signal originates (e.g. 
by detecting Doppler shifts). 



The broad iron line seen from stellar mass black holes ([Miniutti et all 120041) 
provides complementary diagnostics of the strong field region. The profile of 
the line is also broadened and distorted by relativistic effects (Doppler shift 
and beaming, gravitational redshift,...). Combining fast timing and moderate 
resolution spectroscopy will help in testing models and in determining the true 
nature of the underlying variability. As an example, proving the orbital nature 
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Fig. 4. Simulated light curve of a burst oscillation at 300 Hz as seen by the 
RXTE/PCA and XEUS. The burst produces 15000 counts/s in the PCA, and 
400000 count /s in XEUS. The simulations take into account gravitational light 
bending in a Schwarzschild spacetime (one spot and a cosine emission diagram). 
For XEUS, the light curve is computed for two neutron star compactness (M/R=0.1 
and 0.28): the oscillations are directly visible in the light curve even at the highest 
compactness, whereas they are invisible in the PCA data (Barret 2004). 

of the signal could be achieved by measuring independently the frequency, 
from timing, and the radius at which i t is generated from continuum or line 
spectroscopy ( van der Klisl . 



For waveform fitting of X-ray burst and persistent oscillations, a 10 m 2 in- 
strumentation would enable the signal to be detected over one cycle (it will 
also be very sensitive to the harmonic content of the signal). By fitting the 
waveform, it will be possible to investigate the spacetime around the neutron 
star, and simultaneously constrain its mass and radius, and hence constrain 
the equation of st ate of its high density core with unprecedented accuracy 
( Strohmaver . 20041 ) (see Fig 4). 



3.2 XEUS lining up in the ESA Cosmic Vision 



The ESA Cosmic Vision exercise has defined key goals to be addressed in the 
science program for the time period 2015 to 2025. Among those goals, the 
study of matter under extreme conditions (gravity in the strong field limit 
and the equation of state of dense matter) requires a large aperture X-ray ob- 
servatory. In parallel, the XEUS mission concept has recently evolved, thanks 
to the prospects of using lig ht weight Si HPO o ptics, capable of reaching 
a 5 arcsec spatial resolution kmm^il B 

progress is being 

made and continued development will ensure that a mirror area of 10 m 2 can 
be flown with a single launch to the second Earth-Sun Lagrangian point (L2) 
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and autonomous deployment ( Bavdaz et al. . 2003). To address specifically the 



Matter under extreme condition goal of the Cosmic Vision program, a dedi- 
cated timing instrument, capable of coping with the count rates expected from 
the bright sources of the sky is now in the baseline focal plane instrumenta- 
tion. The detector implementation described in Barret (2004) is designed to 
provide an order of magnitude improvement for timing studies over RXTE, 
while providing more than a factor of 5 better spectral resolution, making 
possible to combine timing and spectroscopic measurements. 



4 Conclusions 



RXTE has opened a brand new window on fundamental physics. Gravity in 
the strong field limit and matter at supra-nuclear density, which constitute 
the core of the RXTE legacy, have been identified as key science goals to be 
pursued in the science program of ESA after 2015. XEUS, which is foreseen as 
a potential mission for this time period, now has a dedicated timing detector 
in the focal plane baseline, with capabilities matching the requirements for 
a follow-up instrumentation to RXTE. XEUS has therefore a great potential 
for delivering unique science on strong gravity and dense matter, building on 
what has been learned in the last few years with RXTE. 
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